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ABSTRACT 

Weak lensing is the distortion (polarization) of images of distant objects, such as high 
redshift galaxies, by gravitational fields in the limit where the distortion is small. Grav- 
itational potential fluctuations due to large scale structure cause correlated distortions 
of the images of high redshift galaxies. These distortions are observable with current 
5-H large telescopes and instrumentation. 

In a Friedmann-Robertson- Walker (FRW) metric I calculate the weak lensing pat- 
tern due to large scale structure for arbitrary Slo and zero cosmological constant A. 
For a given cosmological model, specified by Slo and a power spectrum of density fluc- 
tuations, I calculate the statistical properties of the polarization field for an arbitrary 
j redshift source distribution in a simple closed form. 

I It is shown that for low redshift z of the sources, the polarization amplitude is 

proportional to S7o, while at higher redshift the polarization measures the value of 
CO where z is the characteristic source redshift. Moreover, the statistics of the polarization 

field are a direct measure of the power spectrum of density fluctuations. 

ON Key words: Cosmology: theory - gravitational lensing - large scale structure of the 

Universe 
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1 INTRODUCTION 

Gravitational lensing provides a direct probe of mass fluctuations in the Universe. In strong lensing by individual galaxies 
or clusters there is a significant distortion of images even to the point of forming multiple images. This was first proposed 
by Zwicky 1937 and was first observed by Walsh, Carswell & Weyman 1979. Since then, a large number of strong lensing 
systems have been discovered (see f.ex. Schneider, Ehlers & Falco 1992). Weak lensing, on the other hand, is generated by 
sufficiently weak gravitational fields and manifests itself in a distortion of images of distant objects, such as high redshift 
galaxies. A distant galaxy seen behind a high density region will be elongated tangentially relative to the line of sight to the 
cluster center. The same galaxy will be elongated radially relative to the void center if the light rays pass through a void. As 
the light from a galaxy travels through high and low density regions towards the observer the galaxy image will be distorted 
in random directions. 

Weak lensing is therefore a measure of the amplitude and coherence of the cosmological gravitational field and as such is 
a measure of the mass distribution in the universe. The important point about weak lensing is that although on average the 
image distortions will cancel out, within a given small patch of the sky it will induce a preferred direction. 

We parametrize galaxy shapes by the complex eccentricity \ = — eexp(2i0) where e is the eccentricity of the image and <f> 
is the position angle of the major axis. We make the assumption that in this patch of the sky there is intrinsically no preferred 
direction, i.e. < Xtrue >= . The polarization p is a measure of the induced complex eccentricity of an image. If we have an 
image with true complex eccentricity Xtrue, the observed complex eccentricity Xobs will be 



where * denotes complex conjugate. In the presence of weak lensing, the observed mean value of the complex orientation 
< X > is a measure of the polarization p due to weak lensing. < x >~ V- The polarization due to weak lensing is thus an 
observable quantity. 

In an individual cluster the polarization field carries information about the depth and shape of the gravitational potential. 
Weak lensing in clusters has been detected in a large number of systems, e.g. Tyson, Valdes & Wenk 1990, Kaiser & Squires 
1993, Seitz & Schneider I994,Smail et al. 1994. 

If we focus on the polarization due to large scale structure rather than by individual objects, it can be shown that 
the statistical properties of this polarization field carries information about cosmological parameters such as the value of the 
cosmological density parameter fio and the power spectrum of density fluctuations P(k). The first suggestion that weak lensing 
could be used to probe the large scale mass distribution was made by Gunn 1967 and the first observational investigation 
by Kristian 1967. Later fairly shallow searches by e.g. Valdes, Tyson & Jarvis 1983 were also unsuccesful. Recent surveys, 
e.g. Mould et al. 1994 (M94), have yielded marginal detections. These are small, very deep pilot surveys. Their importance 
is that they demonstrated that it is possible to control the systematic errors to the point that larger surveys under excellent 
atmospheric conditions can detect a signal. In particular the polarization signal in the Sloan Survey (Gunn & Knapp 1994) 
will offer the opportunity to map in this way the gravitational field of the Local Supercluster, as has recently been suggested 
by Gould & Villumsen (1994). 

Weak lensing calculations have been performed by a number of authors, e.g. Blandford et al. 1991 (B91), Miralda-Escude 
1991, Kaiser 1992, M94. For a universe with Qo = 1 and a density contrast of approximately unity on a scale of 800kms -1 , 
they predict a mean polarization of a few percent in a deep field of a few square arcmins. This is observable in a modest 
observational program on a large telescope under excellent atmospheric conditions. 

In this paper I calculate the basic statistical properties of the polarization field in a universe with an arbitrary value of 
fio, assuming A = 0. Instrumental in this calculation is the use of the comoving angular diameter distance as the fundamental 
variable, which makes it possible to extend previous calculations for Einstein-de Sitter universes in deriving closed form 
expressions for the various relevant quantities. The calculation is for an arbitrary source redshift distribution, whereby it is 
assumed that this distribution is the same all over the sky. To be able to perform these calculations I have to introduce a 
few technical approximations. Firstly, I use the Born-approximation, which says that the distortion of high redshift galaxy 
images can be calculated as an integral along the geodesic in a universe with no density fluctuations. In addition, it is also 
necessary to assume the plane-parallel approximation. The latter involves the assumption that the only waves that contribute 
significantly to the lensing are the ones nearly perpendicular to the line-of-sight, if the density field is decomposed into plane 
waves. Finally, also for fio 1 Universes I assume I can carry out a Fourier analysis. Such an analysis is strictly valid only 
for wavelengths much smaller than the radius of curvature of the universe, and therefore contains the implicit assumption 
that the contribution of horizon-scale perturbations is negligible. In practice this puts constraints on the power spectrum 
behaviour at low wave numbers. In general, a power spectrum that runs like P(k) cx k n , with n > 0, at small k will behave 
well. Formally the polarization will be finite for n > —2 but the assumptions in this paper will no longer be valid on large 
scales. 

Some of the present paper has been covered previously by B9I , Miralda-Escude 1991, Kaiser 1992, M94. However, they 
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restrict themselves to the simpler problem of an Qo = 1 Universe. I show that the most important cosmological effect is that at 
low redshift the polarization is nearly proportional to Qo- At higher redshift the polarization amplitude is nearly proportional 
to the value of the cosmological density parameter at the characteristic redshift of the sources. Polarization measurements are 
therefore direct measurements of Qo- Moreover, the dependence of the polarization on the redshift z might contain interesting 
information on A. The equivalent calculation for auniverse with a non-zero cosmological constant will be dealt with in a future 
paper. 

An introduction to weak lensing by large scale structure is given in §f . This if followed in §2 by the calculation of the 
polarization field on the sky. The statistical properties of the polarization field are derived in §3. A summary of the equations 
for the polarization field and the basic statististical properties is given in §4, while in §5 I finish with a short discussion of the 
results. 



2 CALCULATING THE POLARIZATION FIELD 



2.1 Basics 



I want to calculate the evolution of the the magnification and shear (distortion of images of high redshift galaxies) in an FEW 
Universe which is not necessarily flat, but that has A = 0. I use the Friedmann equations for the evolution of the expansion 
factor a, the cosmological density parameter fi(a), the Hubble parameter H(a), the comoving angular diameter distance x and 
the affine parameter A. The normalisation is: a = 1 is the present time, t = to, Ho = H(a = 1) = 1, fi(a = 1) = Qo- The speed 
of light is c = 1. The absolute value of the Hubble constant does not explicitly enter into the results. All the equations are 
written with x, the comoving angular diameter distance as the independent variable, x will henceforth be referred to simply 
as "distance". 

The derivations are an extension and generalisation of B91 and M94 for a weakly perturbed FEW Universe. It is not 
necessary to make simplifying assumptions beyond those presented for aflat universe. We define a distance polar D = (Di, D 2 ), 
which is a generalised proper distance. Both components of D are in general complex numbers. The real part of D\ represents 
the distance while the imaginary part measures the image rotation. As shown in B91 the image rotation is negligible and D\ 
will be treated as a real number. D 2 measures the complex rate of shear. In a homogenous universe D\ is the angular diameter 
distance Dos, while D 2 = 0. We can write the magnification AM and polarization p in terms of D . 

These quantities are first order approximations. 

The fundamental equations for the distance polar D are given by B91 and by Seitz, Schneider, & Ehlers 1994. 

d 2 6D 1 



d\ 2 
d 2 D 2 



HhD 1 = bUDx (3) 

dx2 -KD 2 = FDl (4) 

Here 1Z and.7 7 are derived from the Ricci tensor and the Weyl tensor. In this approximation they can be represented by second 
derivatives in comoving coordinates of the comoving Newtonian gravitational potential $. 

Here are coordinates transverse to the direction of the ray. In the plane parallel approximation we can add the second 

derivative of the potential along the ray. 

ft = -a" 5 V 2 $ (7) 

The quantity 81Z is the perturbation of 1Z relative to a homogeneous universe. Because there is no gravitational shear in a 
homogeneous univers, 6T = T . In these equations D\ will be real, so the image rotation will be zero. We set D\ in Eqs.( 3, 
4) to be the unperturbed distance so 

D 1= ax ; SD t = aSx t ; D 2 = aSx 2 ; AM = -2^- ; p = 2^ (8) 
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Here (8x1, 8x2) is the longitudinal and transverse displacement of the beam relative to the beam position in an unperturbed 
universe. Notice that 8x1 is a complex coordinate. We can then rewrite the fundamental equations in terms of x, the coordinate 
along the unperturbed ray. 



d (a8x 
dA 2 " 

d 2 (0,8x2 
d\ 2 

— 2 



— — TZaSxi = 81Zax (9) 
1Za8x2 = Tax (10) 

(11) 



3 QH 2 3 Q 



2 a 2 2 a 

811 = -^Q a~ 5 %) (12) 

9 2 $ 9 2 $ 9 2 $ \ 3 _ 5t .„ 

df 2 d v 2 d£d v I ~ 2 



= -a" 5 ( — " -JTT + 2i— ) = --Qoa- 5 ^'(x) (13) 



8(x_) and .T 7 ' are the density contrast and the shear in the comoving gravitational potential at comoving position x_ at the 
appropriate epoch a = a(x). The source of the magnification is the density fluctuations 8(x_) while the source of the shear of 
the images is the shear T in the comoving potential. 

T is written in this way for notational consistency with 81Z. For a flat universe both 8(x_) and T'{x) are proportional to 
the expansion factor a in linear theory. A set of useful relations for the calculation follows. 

A(t) = a(t')dt' =* A(a) = f (14) 



H(a' 

*(*) = ^ *(«) = f 777^77777 (15) 



«(*') Ja H ( a ') a 

dx —1 d\ —1 dx 1 d 2 x H(a) 

~da~ = H(a)a 2 ' ~da~ = H(a) ' Ix = "a 2 ' rfA 2 = a 3 ^ ^ 

J ff(a) = a- 1 (l + Q (l/a-l)) 1/2 ; = Q + (1 - «o)a (17) 

g = f + (l-Qo)a; Q(a)I7 2 (a) = f (18) 



2.2 Calculation for Single Source 



Let me rephrase our standard equation for a source at distance x. I will only show the calculation of 8x1 since the calculation 
for 8x2 is nearly identical. 

We first change the dependent variable from the affine parameter A to the distance x. The we reduce as much as possible 
the occurences of the expansion factor a and functions of a. 

d 2 (a8x\) ( dx\ 2 d(a8x\)d 2 x 3^ _ 4 . 3^ -i r , n 

-^{d\) + ^dx- A dV + 2°° a ^="2 fioa =* 
c^a^) a _ 4 + 2 Ma6x 1 ) a . 3 + 3 Qoa - 4fei = _£ Qoa -4 %)x ^ 
ax 2 ax 2 2 

d 2 (a8x\) n d(a8xi) „. . 3„ r 3„ r . . 

dx 2 dx 2 2 

a ^ - 2i7(a)a — + It + (1 - Qo) a J 8x1 + 

2 f — _ff (a)a 2 8x\ + a 1 ^ aH(a) -\ — Q,o8x\ = Qq6(x)x =>• 

V. / 2 2 

d 2 8x\ .„ r 3„ <5(x) 

ax 2 2 a 

Equivalently the equation for the shear is 

d 2 8x2 ,~ x r 3^ J 7 ' (x) , . 

-— - + Q - 1 )«z 2 = --«o— — x 20 
dx 2 2 a 



Weak Lensing 5 



The boundary conditions are 

fin =0 ; ^P- = ;fix 2 =0 ; ^ = ; for x = (21) 

(IX (IX 

We have not been completely successful in getting rid of a, we are left with one power of a. It is, however, not profitable 
to express a in terms of x and fio in these equations. We also see that the solutions to the homogeneous equations change 
character depending on whether the universe if open, flat, or closed. The solutions are exponential, linear, or trigonometric 
functions depending on the value of fio- 

With these boundary conditions the general solutions are 

2 fix! r %o *'(* - *') sin [("° - !) 1/2 (* - x ')] , > f o ^ 1 

= -3O / , A 777 dx Ior "o > 1 

x Jo a \ x ) x {Qo - 1) ' {x - x') 

f x 6(x') x' (x — x') , , „ 
= - 3 / -7=4 — " ~ dx for Q o = 1 

Jo <H X ) x 

n ( x fi(x') x'(x - x') sinh [(1 - Q ) 1/2 (x - x')] 
= -3Q / -7=77 — 11 tp, -dx' for Q < 1 22 

Jo < x ) x (l-Q ) 1/2 (x-x') 

2fix 2 f ^'(x') x'(x - x') sin [(Op - I) 1 ' 2 (x - x')] , 

= -3O / — t~ 7\ 777 dx Ior "o > 1 

x Jo a \ x ) x («o - 1) 1 (x - x<) 

= -3 f X ' ( X ~ X 'h x' forQ = l 

Jo a ( x '> x 

[ x IF'(x') x'(x - x') sinh [(1 - fi ) 1/2 (x - x')l , 

= -3Q / \> ( ' l - tjt^ -d* for «o < 1 (23) 

Jo < x ) x (l-Q ) 1/2 (x-x') 

This can be written in a more general way 

— AM = "3"o^ Jjj^x' (1-7) Jo [(fio-l) 1/2 (x-x')] dx' (24) 

V = J o ^) x ' J ° [("o-l) 1/2 (x-x')] dx' (25) 

This is the general solution for a single source at distance x. Here jo(y) = sin («/)/«/ and Jo(0) = 1. Remember that jo(iy) = 
smh(y)/y. 



2.3 Calculation for a Distribution of Sources 



In a normal observational situation, the sources are not all at the same distance. We take a distribution of sources ra(x) so 
that 

2/n 

dx n(x) = 1 (26) 


The upper bound is given by the horizon distance. A potential fluctuation at distance x' will affect only sources at larger 
distance x > x'. The magnification and polarization measured at a given position on the sky is the mean of the polarizations 
of all the sources at different distances at that position. So a density perturbation at distance x induces a polarization for all 
more distant sources at the same position in the sky. Thus the observed magnification AM = —28xi/x and p = 8x2 /x are 

-3Q / ' " dx'^-x'w (x',Q ) (27) 
Jo < x ) 
l>2/n ,pf/ /\ 

-3Q / dx'^-4^-x'w (x',Q ) (28) 
Jo a ( x ) 



-AM 



(x') = I dx n(x) ( 1 — — ) — "4- — — for fi > 1 

K ' ' V x J (fio -l) 1/2 (x-x') 

2 / x' 
dx n(x) (l ) for fio = 1 



((Qo-l) 1/2 (x-x')) 
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2/n ° ( x'\ sinh((l-Q ) 1/2 (x-x')) 

dx n(x) 1 - — 777 for Q < 1 (29) 

or written in a more general fashion 

/•2/n / ;\ 

w (V, Q ) = / dx n(x) 1 - — Jo ((Oo - if' 2 {x - x')) (30) 



This is the general solution for an arbitrary distribution of sources. This is useful if we know, or assume, a mass distribution 
of scatterers. It is also the basis for a statistical distribution of scatterers, ie. if we assume a cosmological model. 



3 CALCULATING STATISTICAL PROPERTIES 



3.1 Polarization Correlation Function 



I now look at the statistical properties of the polarization field on the sky and I do the general case with a distribution of sources. 
The situation where all the sources are at the same distance is just a special case. First I calculate the two-point correlation 
function of the polarization and the magnification. I do a plane wave decomposition of the density field in the universe. The 
Fourier transform of the density field and of T' differ only by a phase factor, and in the plane-parallel approximation the 
two-point correlation functions are identical (B91 and M94). The plane wave decomposition in curved space is only meaningful 
for wavelengths much smaller than the radius of curvature. We also assume the plane-parallel approximation, i.e. that only 
waves nearly perpendicular to the beam are important. Then we can write 

d 3 k ~ 

7 — 77" ^ (k.) exp(zfcx) = 8(x_) 
( 27r ) 

d 3 k ~ &2 
7 — ^8(k) exp(2iip) exp(ikx_) = T (x_) ; tan(^>) = — (31) 
(27r) 3 k\ 

In equations 27 and 28 we see that the expansion factor a enters explicitly. The terms 8(x_')/a and T'{x^~)ja are evaluated 
at the epoch at which the lightray is at that position. Both quantities are expected to evolve in the same way with x since 
they are both second derivatives of the potential. In a flat universe where the growth of density fluctuations is well described 
by linear theory these terms are equal to 8o(x_') and Tq(x!) which are just the values at the current epoch, i.e. a = 1. For a 
non-flat universe the situation is, as usual, not so simple. We set 

^ = f(x(a), fi ) x So(x) ; EE f(x(a), Q ) x To{x) (32) 

In general, / will be a function of position but we are going to assume that the large scale structure evolves simply by 
increasing the density contrast so / = /(a,fio). Notice that this parametrization does not require linear growth of density 
fluctuations. The quantity / will be discussed in more detail in §3.5 . 

We first calculate the two-point correlation of the polarization C pp (9). This is the same as the magnification two-point 
correlation function. The Fourier Transforms of the magnification and polarization differ only by a phase factor. The amplitudes 
are the same and thus the power spectra and two point correlation functions will be identical. Two direction vectors x_ and 
are separated by an angle 9 <C 1, x_ = x_x, and x_' = x_'x' . C vv (9) can then be calculated as the average over the sky of the 
product of the polarizations measured in directions separated by distance 9. For simplicity we let one direction vector be the 
z-axis. 

C pp {9) = (AM(i)AM(i')) = (p(i)p*(i')) 

1o / dx f(x) xw(x)8o(x_) / dx' f'(x') x'w(x')8o(x_') 
Jo Jo 

2 f 2/n ° f d 3 k ~ 

9Q / dx f(x) x w(x) I - — tt8q (k) exvdkx) x 

Jo J ( 27r ) 

2/n ° i i i i if d 3 k ~, , , 
dx f (x ) x w(x ) / - — 77^o(A ) exp( — ikx ) 

o J ( 27r ) 

= 9Qq / dx f(x) x w(x) I dx' f'(x') x' w(x') I d 3 kP(k) exy>(i(kx — kx')) (33) 
Jo Jo J 

We have used the definition of P(k), which is evaluated at the present epoch, so that 
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(«o(*)*o*(*')> = (27r) 6 P(k)8(k-k') (34) 

We now need to use the small angle approximation 9 <C 1 and the plane-parallel approximation, that only waves nearly 
perpendicular to the line-of-sight contribute to the polarization. This enters in the fi integral where the limits of integration 
are extended to ±00. I set 

k_ = A;(cos cf>, sin cf>, fi) ; x_ = x(0, 0, 1) ; x_ = x' (9, 0,1) ; x_ — x_ = ( — x'9, 0, x — x') (35) 
d 3 kP(k) exy(i(kx — kx')) 

dkk 2 P(k) I d(f> exp( — ik cos (f>x'ff) / dfi exp(ikfi(x — x')) 
Jo J — co 

f CO 

= 6(x-x')4:Tr 2 dkkP(k)J (kx9) (36) 
Jo 

We insert this result into the equation for cpp. 

/•2/rio /"co 

C pp {9) = 36Qot 2 / dx f 2 (x)x 2 w 2 (x) / dkkP(k)Jo(kx9) (37) 
Jo Jo 

We see that the two-point correlation function is a first moment of the power spectrum of density fluctuations. For a source 
at distance x, only waves with wavenumber less than (x9)~ 1 contribute significantly to the correlation function. 

We wish to relate the two-point polarization correlation function to the two-point correlation function of the mass {(r). 
For this we can use the relationship between P(k) and {(r), they are Fourier Transforms of each other. 

p W = / dyy 2 t(y)Mky) (38) 

/7r Jo 

When we insert this relationship we get for C pp (9) that 

/•2/rio /"CO /"CO 

Cpp{9) = 18Qo / dx f 2 (x)x 2 w 2 (x) I dyy 2 {(«/) / dkkJo(kx9)jo(ky) 
Jo Jo Jo 

2/n yco 

dx f 2 (x)x 2 w 2 (x) I dyyS,{y) (y 2 - (x9) 2 ) 

Jx6 

poo 

dx f 2 (x)x 2 w 2 (x) dy£ (y / (x9) 2 + yA (39) 

J ^ ' 

= 9Qo / dx f(x)x 2 w 2 (x)W(x9) (40) 
Jo 

For a source of distance x, the contribution to C pp (9) can be written as an integral over the correlation function on scales 
larger than x9. The quantity W(x9) is the two-point correlation function of the surface density distribution evaluated at angle 
9. As a simple example of what might be observed is seen in Figure 1 I show [C pp (9)] 1/2 for a CDM and an HDM power 
spectrum, Bardeen et al. 1986, where the powerspectra are normalized so that the rms density fluctuation in a sphere of radius 
8/t -1 Mpc for h=l/2, equals unity. For simplicity it is assumed that all sources are at redshift 2=3/4 and fio = 1. 

The lower curves (CDM solid, HDM dotted) are the curves for the correlation function on a scale up to 20 arcmins. An 
angular scale of 10 arcmins corresponds to a linear scale of 4.2 ft -1 Mpc at that redshift. On that angular scale both CDM 
and HDM have a correlation amplitude of (0.02) 2 . On smaller scales CDM greatly increases that amplitude while for HDM 
the amplitude stays constant. 



3.2 Polarization Power Spectrum 



Let us look at the power spectrum of polarisation fluctuations Q(zu). The power spectrum of polarization fluctuations is the 
Fourier transform of the two-point correlation function C pp (9). Notice that this definition of Q is different from the definition 
in M94 and Gould & Villumsen 1994. In those papers Q = Q(k) where k is an inverse length while in this paper Q = Q(zu), 
zu is an inverse angle. If the all the sources are at distance x the two definitions differ by a factor x 2 . 

i Z" 00 

Q(zu) = — d9 9C pp (9)Jo(zu9) 
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Figure 1. Square root of polarization correlation function C pp (6) and rms polarization fluctuations for CDM (solid lines), and HDM 
(dotted lines). All sources at z = 3/4, Qq = 1. Upper curves are rms fluctuations, lower curves are the square root of correlation 
amplitude. 



18ttQo 



18ttQo 



18ttQo 



2/n 



dx f 2 (x)w 2 (x)x^ 



2/n 



dx f 2 (x)w 2 (x)x^ 



dkkP(k) d0 0J o (kx0)J o (zu0) 



dkkP(k) 



8(k — zu I x) 

x 2 k 



2/n 



dx f 2 (x)w 2 (x)P(zu / x) 



(41) 



Notice that in this derivation I have used 8 as symbol for the Dirac delta function. We can use the power spectrum of potential 
fluctuations P#(fc) which is related to P(k) 



V 2 $(x) 
Q(zu) 



18ttQo 
8ttG 2 



$(i) 



4ttG~ 



it) 



l 67r 2 G 2^) 



dkf 2 (zu/k)w 2 (zu/k)P(k)k- 



dkf 2 (zu/k)w 2 (zu/k)Pz(k)k 2 



(42) 
(43) 

(44) 



This equation says that the contribution to the power spectrum of polarization fluctuations at inverse angle zu from 
structure at a distance x comes from the power spectrum of density fluctuations at wavenumbers larger than zu/x. 

In Figure 2 I show the standard power spectra of density fluctuations P(k)k 2 for CDM and HDM, Bardeen et a,l. 1986 
(Solid lines). The normalisations of amplitude are arbitrary. On the same plot is shown Q(zu)zu for all sources at distance 
x = l,ie. z = 3 (Dotted lines). This is done in order to directly compare the powerspectra for density fluctuations and for 
polarization fluctuations. In this case zu = k. I have multiplied P(k) by two powers of the wavenumber and Q(k = zu) by one 
power to demonstrate on what scales the density fluctuations and polarization fluctuations are generated. We see that the 
polarization fluctuations are generated at significantly larger scales than the density fluctuations. The range of scales that 
contribute to the polarization is also larger. 

We can get some insight into equations 43 and 44 by considering a simple calculation for a flat universe. Suppose all the 
sources are at distance xo and P(k) can be locally approximated by a power law of slope n then we have approximately that 
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Figure 2. Power spectra of density fluctuations P(k)k 2 , and polarization fluctuations Q{m = k)k for CDM and HDM for sources at 
distance x = 1. All absolute amplitudes are arbitrary. 



^oc r dk{\-^-\ (45) 
This integrand peaks at 

— = k p = — =>• x p = x (46) 



2 — n 7 xo 

This says that the peak contribution to the polarization occurs at distance x p from the observer. Let us consider CDM as an 
example. If we are looking on the Harrison-Zeldovich tail,ra = 1, then x p /xo = 1/3 meaning scales closer to the observer than 
to the source are most effective. If we look at galaxy scales, n ~ —2, then x p /xo = 2/3, most of the power is generated nearer 
the source. For positive slope of P(k), more of the power is on smaller scale, which pushes the most effective lens distance 
towards the observer. For a negative slope of P(k), the opposite is true. 

How meaningful is the simple / parametrization? We see from Eq. 44 that the power spectrum Q can be derived as 
an integral over the power spectrum of potential fluctuations with a high pass filter w. Potential fluctuations are weighted 
towards large scales so the simple parametrization is expected to work well. Notice again, that this parametrization does not 
require linear theory. 



3.3 Polarization Variance 



Another important statistical quantity is the polarization variance on a scale 9 using an axially symmetric window function 
W 1 (y/8), where 

p CO 

2tt / dyyW 1 (y/0) = 1 (47) 
./o 

W 1 (w6)= [ ^LyW 1 (y/6)Jo(yw) (48) 

/ Z7T 
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The variance of a homogeneous field on a scale 9 is the mean square amplitude of the field when convolved with a filter 
function Wi(y/9). This variance can be evaluated simply as 



<Tp(6) = 2tt / dwwQ(w)Wf (w6) 



o 

= 36Qo*- 2 / dx f 2 (x)x 2 w 2 (x) dkkP(k)W 2 (kx9) (49) 
Jo Jo 

where we have used Eq. 41 and that k ■ x = zu. For the special cases of a tophat filter and a Gaussian filter of radius 9 we get 
a 2 p (9) = 36Q^7r 2 ^ ' ° dx f(x)x 2 w 2 (x) dkkP(k) ^^fl^ ) ( 50 ) 

°pW = 36Qo7r 2 / dx f 2 (x)x 2 w 2 (x) I dkkP(k) exp( — k 2 x 2 9 2 ) (51) 
Jo Jo 

The inner integral in Eqs. 49, 50, and 51 is just w/2 times the mean square density fluctuations cr 2 , on a sheet perpendicular 
to the line of sight with the appropriate window function. Notice that <r M is the volume density fluctuations evaluated on a 
sheet, not the surface density fluctuations. Thus we get that 
/•2/fi 

(T 2 p {9) = 18fio7r 3 / dx f 2 (x)x 2 w 2 (x)a 2 l (x9) (52) 
Jo 

Neither the two-point correlation function C pp (9) nor the variance cr 2 , depend on the statistics of the polarization field. Figure 
1 also shows a p for CDM and HDM for the same parameters as in §3.1(upper curves). The filter function is a tophat. On a 
scale of 5 arcmin, which is the scale measured by M94, the CDM prediction is an rms amplitude of 3% while HDM predicts 
2% for Qo = 1. 



3.4 Gaussian Polarization Statistics 



If the field of density fluctuations is a Gaussian random field then the polarization field will also be a Gaussian random field. 
Polarization results from scattering of lightrays off the density fluctuations. If the density field is not Gaussian then if the 
scatterers are uncorrelated, by the Central Limit Theorem, the polarization field will become Gaussian for a sufficiently large 
pathlength. We know that the density field is correlated but if the pathlength is much larger than the correlation length of 
the density field the Central Limit Theorem will still imply that the polarization field will be Gaussian. 

In that case Q(zu) fully specifies the statistics of the polarization field. Eq. 41, or equivalently Eqs. 43, 44, is the 
fundamental equation for the statistical study of weak lensing. Given a known source distribution, then from the power 
spectrum of density fluctuations and Qo we can predict the polarization field. The fact that the polarization field is Gaussian 
means there will be specific predictions for the higher order correlations of the field given the two-point correlation. The 
distribution of polarizations will be Gaussian. If we have observed a polarization field, i.e. < x >, as a function of position 
on the sky, then if this field does not have Gaussian statistics or the higher order correlations do not agree with predictions 
based on the two-point correlation function, then we can falsify the weak lensing hypothesis for this < \ > field. There is a 
caveat to this statement, single scattering events from individual clusters can be stronger than the polarization in the field. 
Thus in the tail of the distribution we can be dominated by single scattering events and the polarization distribution need 
then not be Gaussian. 

If the polarization field is a complex Gaussian field, then the amplitude field will have a Rayleigh distribution. In that 
case, the most likely measured amplitude will be the rms amplitude. Thus for M94 the most likely polarization amplitude in 
an unbiased, flat CDM universe will be 3%. 



3.5 Estimate of Curvature Effects 



Let us look at the effects of varying Qo. For simplicity we put all sources at distance x, or equivalently redshift z. This makes 
the weightfunction w simple 



(x', Q ) = ( 1 - - ) jo ((Qo - l) 1/2 (x - x')) (53) 
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We look at the rms polarization <r p (0) as a function of fio for a given source redshift z. At low redshift, z <C 1, we have x = z 
for any value of fio and the weight function w is independent of fio- That means <r p (0) is proportional to fio for a fixed source 
redshift. This proportionality makes physical sense since then the polarization is proportional to the potential well depth. If 
we go to higher values of z three effects enter that all pull in the same direction. For a given z, as we lower fio, the distance 
increases, the weight function w increases, and / increases. The distance is larger because normalised to the present epoch 
the universe is older and therefore the lightray has travelled further. Secondly, w has increased because x has increased and 
the lightrays are more divergent (for fio < 1) or less convergent (for fio > 1). The different cosmological models also enter in 
the terms 8(x_')/a and T'{x^)ja evaluated at the epoch at which the lightray is at that position. Both quantities are expected 
to evolve in the same way with x since they are both second derivatives of the potential. In a flat universe where the growth 
of density fluctuations on relevant scales is well described by linear theory these terms are just equal to 8o(x_') and Tq(x^) 
which are just the values at the current epoch, i.e. a = 1. For a non-flat universe the equations for / are more complicated, 
however the validity of the linear theory estimate of / is not affected. In linear theory (Peebles 1980) the amplitude of the 
growing mode, and therefore /, can be calculated for any value of fio- We can get a rough estimate of / from linear theory 
where the logarithmic growth rate of the growing mode is approximately Qq' 6 . Thus we get 

/ « a ("o 6 -i) ~ (i _ z) « - 6 -i) «(i+[i_ n° 6 ] z) for z < I (54) 

We see from this simple parametrisation that / decreases with fio and increases rapidly with distance. This effect is the most 
important of the three effects. 

B9I indicate from N-body simulation for Cold Dark Matter and Hot Dark Matter that for the purposes of weak lensing 
calculations it is a good approximation to use linear theory for the growth of density fluctuations. There is, however, a finite 
probability that light rays will pass close to a dense cluster where the deflections will be strong. In this case the formalism 
will break down and we need to talk about strong lensing. However, the fraction of sky covered with strong lenses is small 
and this problem will be ignored in the rest of the paper. 

From Eq 52 we see that the amplitude of the curvature effects is independent of P(k). We can thus show simply the 
curvature correction C, which we define as the rms polarization fluctuations at zero separation as a function of fio and the 
redshift of the sources z measured relative to the rms fluctuations for fio = 1. Here we have for illustrative purposes assumed 
that all the source are at the same redshift. 

C(ii °- Z)= MO; Z ;iio=l) (55) 
The curvature correction can be interpreted as the value of the density parameter that we would infer from observations if we 
naively assumed that the universe were flat (Qo = I), and that the growth of density fluctuations were linear in the expansion 
factor a. 

In Figure 3 we show the curvature correction C(Qo, z) as a function of z for fixed values of fio ranging from a very open 
universe (Qo = 0.2) to a very closed universe (Qo = 2.0). At low redshifts we see that C(Qo,2 <C 1) ~ Oo but at higher 
redshift there are significant extra curvature effects. For open universes the curvature correction is an increasing function of 
source redshift z while C decreases with z for a closed universe. In other words, lensing amplitudes will approach the values 
for a flat universe. A simple way of looking at this is to see that as we go to higher redshifts, Q(z,Qo) will approach unity. 
This is demonstrated by the dotted curves where we have plotted Q(Qo,2) on top of the curvature correction for the same 
values of fio in Figure 3. For a given value of fio these curves lie nearly on top of each other showing that what we are really 
measuring is the value of Q at the epoch when the photons were emitted from the source. In other words 

c(n ,*)«n(no,*)= "° + " z (56) 

I + \Iq z 

That means if we know the source redshift and the value of Q at that epoch we can infer the current value fio from the 
equation. 

Qo = i + (i n -n ( «))« < 57 > 

We can thus write an approximate equation for the rms polarization fluctuations in terms of source redshift z. 

3 / 1 _|_ z \ 2 f°° 

a 2 J0) = C pp (0) x 36tt 2 (2 [f - (I + z)~ 1/2 } ) QjU / dkkP(k) for z < 1 (58) 

VI + IIqzJ J 

f CO 

<Tp(0) = Cpp(0) « 36ttV(I + 2(1 - tt )z)tt 2 o / dkkP(k) for z < I (59) 

Jo 

In the same way we can write C pp (9) and cr p (d) in terms of the source redshift z for a given value of fio- Since x(Qo, z) is a 
decreasing function of fio for fixed z we find that as we lower fio, both C pp (9) and cr p (d) will be steeper functions of 9. 

C pp {6,z,Q ) x, tit I" 1 + * Y C pp (8 ■x(Q ,z)/x(n = l,z),z,O = I), z<l (60) 
VI + \loz J 
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Figure 3. Curvature corrections for rms polarization fluctuations. Solid lines shows C(Qo,z) a t different source redshifts for a set of 
values of Qq. The dotted lines show the value of Q at the source redshift. 

<rl{e,z,U ) « Qo (-T-TT-) 2 "It 6 ■ x(Q ,z)/x(Q = \,z),z,U = 1), z<l (61) 
VI + iloz / 

The weak lensing prediction for a zero cosmological constant universe is that observations of the polarization amplitude as a 
function of source redshift will follow a curve like those shown in Figure 3. If the amplitude does not follow such a curve, it 
is an indication for non-zero A. 



4 SUMMARY OF RESULTS 



In a universe with an FRW metric, defined by the current cosmological density parameter fio, and a zero cosmological constant, 
we have a density field 8o(x_), defined at the current epoch with a power spectrum P(k). It is assumed that the growth of 
density fluctuations can be simply described by a universal function /(x,Qo)- It is not necessary to assume linear growth 
factors. The shear in the gravitational field is desribed by the terms T' defined in equation 6. The distribution of sources in 
comoving angular diameter distance x is n(x). Other terms are the expansion a, the inverse angle zu, the power spectrum of 
potential fluctuation _P#, and the variance of the volume density on a sheet, W. 

(62) 



/ = 


8(x) 


_ T'(x) 


8o ( x_) a 


T' a (x_) a 


w (x' , Qo) = 


J X 1 


dx n(x) ^1 


—AM = 


-3Q 


/ dx' 8 \ 
'o a ( 







^Jjo U<h - 1 (>-.:■)) (63) 
l-x'w (x',Q ) (64) 



y2/n T'(x') 

-3Qo / dx — ~ x w (x , fio) (65) 
' a(x') ' 

C pp (8) = 36QqT 2 / dx f(x)x 2 w 2 (x) / dkkP(k)J o (kx0) (66) 
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/•2/fi 

= 9Qq / dx f 2 (x)x 2 w 2 (x)W(xff) (67) 
Jo 

Q(w) = 18ttQo / dx f(x)w 2 (x)P(w/x) (68) 

no2 /* co 

= 77? / dkf(wlk)w 2 (wlk)Ps(k)k 2 (69) 
1*2 /r^Q /• co 

<r 2 (6>) = 36Qo?r 2 / dx f 2 (x)x 2 w 2 (x) dkkP(k)W 2 (kx0) (70) 
Jo Jo 

= 36Q^tt 2 ^ 7 ° dx f(x)x 2 w 2 (x) dkkP(k) ( y 2Jl ^^ Tophat (71) 

= 36Qo 71 ' 2 / rfx / 2 (x)x 2 w 2 (x) / dkkP(k) exp( — k 2 x 2 9 2 ) Gaussian (72) 
Jo Jo 

If all the sources are the same redshift z we can write the results in terms of redshift, where C is the curvature correction. 



<t 2 (0) = C pp (0) « 36tt 2 (2 [l - (1 + z)- 1/2 ] ) 3 O 2 , ^ dkkP(k) for 2 ^ 1 (73) 

/» CO 

<r 2 (0) = Cpp(0) « 36ttV(1 + 2(1 - tt )z)tt 2 / dkkP(k) for 2 < 1 (74) 

Jo 

C PP (0,2,Q O ) « "of / + * Y C pp (0 -x(Q ,z)/x(Q = \,z),z,U = 1), z< 1 (75) 
VI + s2o2/ 

<t 2 (0,2,Q o ) « "o 1 + Z V o>(# -x(O ,g)/x(O = l,z),z,n = 1), z£l (76) 
VI + iloz J 

C(fio^) « O(0o,^)= "°, + " * (77) 

1 + I Iq z 



5 DISCUSSION 



The polarization field has been calculated rigorously, in the Born- approximation, for any value of the cosmological density 
parameter fio, but with zero cosmological constant. For a known mass distribution we calculate the polarization for a distri- 
bution of sources in real space. We show rigorously that at low redshift the polarization is proportional to fio- Thus weak 
lensing observations, such as in M94, is a direct measure of fio- This is in retrospect not surprising, since in that situation the 
scales of the universe that we are probing are much smaller than the radius of curvature of the universe. At higher redshifts 
explicit curvature effects manifest themselves in the calculation of w(x) but only weakly. The main curvature effects are that 
the comoving angular distance x is no longer equal to the redshift z and that the growth of density fluctuations depends on 
fio- This formalism permits the comparison of weak lensing observations of large scale structure with theoretical models of 
the universe in an FRW metric. 

In summary, the current formalism permits predictions of the polarization field and its statistical properties for a well 
defined cosmological model in an FRW universe with A = 0. 
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